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^— ^ We present a microscopic theory of single-electron tunneling through metallic nanoparticles con- 

I nected to the electrodes through molecular bridges. It combines the theory of electron transport 

through molecular junctions with the description of the charging dynamics on the nanoparticles. 
We apply the theory to study single-electron tunneling through a gold nanoparticle connected to 
the gold electrodes through two representative benzene-based molecules. We calculate the back- 
ground charge on the nanoparticle induced by the charge transfer between the nanoparticle and 
qq ■ linker molecules, the capacitance and resistance of molecular junction using a first-principles based 

Non-Equilibrium Green's Function theory. We demonstrate the variety of transport characteristics 
that can be achieved through "engineering" of the metal-molecule interaction. 
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Chemically tailored metallic and semiconducting nanoparticles and their assemblies have become the model system 
! for studying the fundamental physics and chemistry of nanostructured materials^ The advancement of the fabrication 
jjy • process using self-assemblyi or biodirected-assembly technique^ with molecular recognition has made devices based 
on single and assembled nanoparticles attractive candidates in applications including single-electronics, novel biosen- 
sors and nanophotonic devices4*2i&4 For single-electron device applications, the possibility of accurate control over 
particle size and density provides a substantial advantage over conventional granular metal thin films and suggests 
possible room-temperature operationi^ Understanding electrical transport in such metal-molecule composite systems 
is therefore an important problem. 

In the "orthodox" theory of single-electron tunneling through a metallic island separated from the electrodes by 
insulating gap^ the physics of Coulomb blockade can be well understood using simple circuit-level theory, with circuit 
parameters including junction resistance, capacitance and charging energy often obtained by fitting experimental data 
and/or using simple electrostatic considerations. In contrast, for the nanoparticle-based single-electron devices, the 
transfer of single electron is achieved by tunneling through the linker molecule connecting the nanoparticles to the 
electrodes and to each otheri^ The linker molecules can be coupled strongly to the electrodes and nanoparticles 
through appropriate end groups, and the tunneling barrier is induced by the energy mismatch of the molecular levels 
relative to the metal Fermi-level £ For the molecular junction, the metallic screening of the applied electric field 
Qv^ ■ occurs over a distance comparable to the size of the molecule 7 ^. Consequently, the capacitance of the molecular 
junction becomes an electro-chemical quantity rather than a geometrical quantity, whose evaluation requires a self- 
£NJ ■ consistent analysis of the charge/potential response within the molecular junctions^ In addition, charge transfer 
between the linker molecule and the metals will lead to fractional charge on the nanoparticle, 7,8 giving rise to an 
intrinsic background charge (even in the absence of charged impurities and gate voltages), which varies in general 
with the applied bias voltage and may affect significantly the current-voltage characteristics^ Better understanding 
of such molecular-assembled single-electron devices therefore requires combining the single-electron tunneling effect 
with a microscopic description of the electronic processes in the molecular tunnel junctions. 

In this paper we describe a microscopic theory of single-electron tunneling through molecular-assembled metallic 
• nanoparticles by combining a first-principles Non-Equilibrium Green's Function (NEGF) theory of molecular trans- 
poriiS with a real-time perturbation theory of the reduced density matrix describing the charging dynamics of the 
metallic nanoparticlesii*^ In the lowest-order perturbation theory in the reduced conductance of the molecular junc- 
tion a = 2 J 1 Rt (Rt is the junction resistance), the theory formally reduces to the rate equation of the "orthodox" 
theory of Coulomb blockade^ but allows us to take full account of metal-molecule interaction at the atomic scale. 
In particular, the background charge, capacitance and tunneling rate can all be obtained from the first-principles 
theory of the molecular tunnel junction. We apply the theory to single-electron tunneling through gold nanoparticles 
connected to the gold electrodes through two benzene-based molecules-biphenyl dithiolate (BPD) and difluorobenzene 
(FBF) molecules and show the the variety of current /conductance- voltage characteristics that can obtained through 
"engineering" of the metal-molecule interaction. 

A schematic illustration of the double-barrier tunneling junction is shown in Fig.^ The system is described by the 
following Hamiltonian: 

H = H L + H R + H I + V + H Ml + H Ah + H Tl + H T2 . (1) 

where H a — ^ fe e^ a au aka and Hi = eucjjCu describe the noninteracting electrons in the two leads (a = L,R) 
and on the metallic nanoparticle (I) respectively. The indices k and I enumerate the electron states of the leads 
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and the nanoparticle. The Coulomb repulsion on the nanoparticle V is obtained from electrostatic considerations as 
V(N) = Eq(N — n x ) 2 , where Eq is the charging energy. The background charge en x can be separated into an intrinsic 
background charge en x o existing at zero voltage and voltage-dependent polarization charge en x = Qp(Vl, Vr) + en x o- 
The junction capacitances are obtained from C a — ^§\r~i where V a ,C a (a = L,R) are the voltage and capacitance 
of the two tunnel junctions respectively. The charging energy is Ec — e 2 /2Cs and the total capacitance Ce = 
Cl + Cr + Cs is the summation of the two junction capacitances and the self-capacitance Cs of the nanoparticle. 
Here the background charge en x is induced by the charge transfer between the nanoparticle and the linker molecules, 
whose value at zero bias voltage gives the intrinsic background charge en X Q. 

We assume the molecules can be described by an effective single- particle Hamiltonian = J2 m e imb\ m bi m , i = l,2. 
The transfer of single electrons is mediated by tunneling through the molecular bridges = E m); j((mio!'] m flfa + 
tmliol m c 'li e + C.C.), where a denotes L (R) for molecule 1 (2). The operators e ±l * keep track of the change of 
the charge on the metallic nanoparticle by ±e^*iii^ The phase operator <p is the quantum-mechanical conjugate of 
the excess electron number N on the nanoparticle. The charging state N is treated independently of the degrees of 
freedom described by the fermionic field operators c u ,cii, which is a good approximation for metallic nanoparticles. 

The two electrodes as well as the nanoparticle are treated as large equilibrium reservoirs with corresponding Fermi 
distribution fi{E) = f(E — fii),i = L,R,I. Note that the partition of the voltage drop between the two molecular 
junctions is determined by the capacitance ratio Vl/Vr = {pi — ^l)/(^r. — Hi) = Cr/Cl and eV = [Ir — ^l- 

Following the standard procedure in NEGF theory of mesoscopic transport^ the current flowing through the left 
contact is given by 

hit) = ^ = ^-Re[J2t mkL G< Lm (t,t)] 



dt 

mkL 
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J ^Tr[G> i (t,t')S<( i ',t)-C?< i (t,t')S>( i ,t)] (2) 

where the Green's functions are defined in the standard manner^ The self-energy operator S^ 5 "^ describes the inter- 
action with the left electrode £^.^(t, if) = J2k ^mkl(t)9L^k (*> O^nfciOO (ffi i s tne corresponding Green's function 
of the bare left electrode) . Similar definition applies to the self-energy operators of the right electrode (R) and the 
central island (I) . Deriving the equation of motion of the Green's function and decoupling the charging dynamics from 

the single-electron dynamics on the metallic nanoparticle {cl I (t)cu(t')e t ^ t "> e' 1 ^} = {clj(t)cii(t')) (e 1 ^" 1 e' 1 ^} , we 
obtain the current as^ 

I L = - 2ei J ^Tl(E,V) J d W 

[f L (E)(l - fi(E - w))C>(w) + (1 - f L (E))fj(E - w)C«(w)}, (3) 

where C < ^ > \w) are the Fourier transform of the correlation functions C < (t,t') = i(e"^* 'e - *^')), 
C > (t,t') = —i{e~ l ^ t ^ > e l ^ t ^). The "renormalized" transmission function is obtained from Tl(E,V) = 
Tr[r L (E,V)G r Mi (E,V)Ti(E,V)G a Mi (E,V)} and the "renormalized" Green's function is G r ^(E,V) = (E ± 0+ - 

H M \ ~ Tf L {a) {E,V) - J dwZ r j {a) (E - w^C^iw))- 1 where C r{ - a \w) are the Fourier transform of C r ^{t,t') = 
±9(±(t-t'))(C>(t,t')-C<(t,t')). 

Eq. can be recast into a form suitable for perturbation expansion 

I L = -ie J dw[aj;(w)C > (w) + a£(w)C < (w)], 

<*1W = 2^"/ dEf L (E,V)f L (E)(l-fi(E-w)), (4) 

«ZM = 2^/ dEf L (E,V)(l-f L (E))fi(E-w), 

where at is the rate for electron tunneling into (out of) the nanoparticle through the left molecular bridge. A 
similar equation can be written down for the current flowing through the right contact Ir = —ie J dw[a^(w)C > (w) + 

a~j i (w)C < (w)} where the tunneling rate oi^ is evaluated from the transmission coefficient through the right molecular 
junction and the corresponding Fermi distributions. The condition for current conservation II, + Ir = leads to 

dw[a + (w)C > (w)+a-(w)C < (w)]=0, (5) 
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where a + (~^(w) — a^'(w) + a^~'(w). Eqs. J2J and J5J are the central result of this paper which separate the 
description of electron tunneling through the molecular bridges from that of the charging dynamics of the nanoparticle 
through the correlation functions C < ( > \ 

The theory of the correlation functions C <( ^ > - ) has been developed by Sch611er,Sch6n and coworkers in a series of 
publications based on a real-time perturbation theory of the reduced density matrix describing the charging state 
N of the metal island, 2U2k where the perturbation parameter is dimcnsionless conductance of the tunnel junction 
ao = h/2e 2 RT- In particular, the correlation function C < ^ > ^ can be expressed as the superposition of contributions 
due to the charging state transition from n to n + 1 as C <( - > ^(?«) = 'Y2 n C <t ^ > \w,n). The current conservation 
condition JSJ can also be shown to hold for each part C < ^ > - ) (w, n)^> J dw[a + (w)C > (w, n) + a~ {w)C < (w, n)] — 0. 
In the lowest order perturbation theory with respect to ao, C <( - > ^ (w, n) are related to the probability P n of finding 
n excess electrons asi£ C < (w,n) = 2TriP n+ i5(w — A„) and C > (w,n) = —2iiiP n 8(w — A„). Here A„ is the energy 
difference of the adjacent charging states A„ = V(n + 1) — V(n) = Eq [1 + 2(n — n x )] . Substituting the above formulas 
into Eqs. Q and (JSJ), we obtain a set of rate equations formally equivalent to the "orthodox" theory of single-electron 
tunneling 

a+{A n )P n - a-{A n )P n+1 = 0, (6) 

Compared to the "orthodox" theory, the above theory allows us to take full account of the microscopic description 
of the electronic processes in the molecular junction. In particular, considerable simplification can be achieved in the 
"wide-band limit" of the electrodes, i.e., if the surface density of states of the electrodes is approximately constant 
in the energy interval of A„. The "renormalization" of the nanoparticle self-energy and transmission coefficient can 
then be neglected, which is true for common simple and noble metals and for nanoparticle size sufficiently large that 
energy quantization can be neglected^ 

Here we apply the theory to single-electron tunneling through a gold nanoparticle connected to two gold electrodes 
using biphenyl dithiolate (BPD) and difluorobenzene (FBF) molecules. Transport through gold-BPD-gold and gold- 
FBF-gold junctions have been recently studied in detail using a first-principles self-consistent matrix Green's function 
theory of electron transport in molecular junctions, 7,8 from which we obtain the voltage-dependent background charge 
en x induced by the charge transfer from the molecule to the nanoparticle and the transmission coefficient for electron 
tunneling through the metal-molecule-nanoparticle junction^ The charge transfer at zero bias gives the intrinsic 
background charge en X Q, while its derivative with respect to bias voltage gives the capacitance of the molecular 
junction^ The calculated dimensionless conductance, capacitance and the intrinsic background charge for BPD 
(FBF) molecular junction are 0.0183 (0.0023), 0.037 (0.043) aF and -0.623 (-0.805) number of electrons respectively. 
Note the total intrinsic background charge on the nanoparticle is the sum of contributions from the two molecular 
junctions. Despite the shorter length of the FBF molecule, the resistance of the FBF junction is much higher than 
that of the BPD molecule due to the less favorable energy-level lineup relative to the gold Fermi-levelA^ 

Given the transmission coefficient and capacitance of the molecular junction, the rate equation © is solved using 
the standard procedure^ from which we calculate the terminal current as 

I L[R) = 27re^P n [a+ (i j ) (A n ) - a~ (jR) (A„_i)]. (7) 

n 

The average excess charge e(N) = e^2 N NPjy gives the net charge transferred onto the nanoparticle due to the 
discrete electron tunneling across the molecular bridges, while the average charge e((N) — n x ) gives the charging 
configuration which determines the electrostatic energy cost for charging the nanoparticle. 

Figs. |21 and 13 show the calculated results for single-electron tunneling through a 10 nm-diameter gold nanoparticle 
connected to gold electrodes through two BPD or FBF bridges at temperatures of 10(K) and 300(i"T) respectively. 
The self-capacitance of the nanoparticle is evaluated using the electrostatic formula for a conducting sphere in front 
of a conducting plane^ which gives a charging energy of Ec = 88(meV). The intrinsic background charge n X Q on the 
central nanoparticle is closer to half integer in the case of the FBF molecule (—1.61) than that of the BPD molecule 
(—1.25), so the Coulomb gap in the current/conductance- voltage characteristics is smaller. Note at low temperature 
the average excess electron (N) on the nanoparticle at low bias equals the nearest integer to n X Q. Since the background 
charge n x varies slowly with the applied bias voltage within the voltage range studied, both the average charge and 
the average excess charge show similar voltage dependence. The peaks in the conductance-voltage characteristics 
correlate with the change in the average excess electron on the nanoparticle at low temperature. 

Fig. 2| shows the calculated results for single-electron tunneling through a 10 nm-diameter gold nanoparticle 
connected to left gold electrode through BPD molecule and connected to right gold electrode through FBF molecule. 
Since Rfbf / Rbpd ^ 1, varying bias voltage leads to stepwise charging of the central nanoparticle at low temperature. 
Here the intrinsic background charge n X Q w —1.43 and Cfbf/Cbpd > 1, the shape of the Coulomb staircase 
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corresponds to case I as discussed by Hanna and Tinkham 10 and is well understood. The small Coulomb gap is again 
due to the fact that n x0 is close to half integer. 

To conclude, we have presented a microscopic theory of single-electron tunneling through molecular-assembled 
metallic nanoparticles, which combines the theory of transport through a molecular junction with the description 
of the charging dynamics on the nanoparticle. This allows us to take full account of the microscopic description of 
the electronic processes in the metal-molecule-nanoparticle junction as well as the strong Coulomb interaction on the 
nanoparticle. 

This work was supported by the DARPA Molectronics program, the DoD-DURINT program and the NSF Nan- 
otechnology Initiative. 
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FIG. 1: Schematic illustration of a metallic nanoparticle connected to the source and drain electrodes through two molecular 
bridges. 
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FIG. 2: Single-electron tunneling through a 10 nm-diameter gold nanoparticle connected to the source and drain electrodes 
through the BPD molecule. Left figure shows the current- volt age and conductance-voltage characteristics. Right figure shows 
the average charge e(N — n x ) and average excess charge eN on the nanoparticle in unit of electron charge e. 
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FIG. 3: Single-electron tunneling through a 10 nm-diameter gold nanoparticle connected to the source and drain electrodes 
through the FBF molecule as in Fig. [5] 
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FIG. 4: Single-electron tunneling through a 10 nm-diameter gold nanopartcile connected to the source electrode through the 
BPD molecule and the drain electrode through the FBF molecule as in Fig. [5] 



